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ABSTRACT 

We study the capture of galactic dark matter particles in the Solar System produced 
by rotation of Jupiter. It is shown that the capture cross-section is much larger than 
the area of Jupiter orbit being inversely diverging at small particle energy. We show 
that the dynamics of captured particles is chaotic and is well described by a simple 
symplectic dark map. This dark map description allows to simulate the scattering 
and dynamics of 10 dark matter particles during the life time of the Solar System 
and to determine dark matter density profile as a function of distance from the Sun. 
The mass of captured dark matter in the radius of Neptune orbit is estimated to be 
2 • 10 15 g. The radial density of captured dark matter is found to be approximately 
constant behind Jupiter orbit being similar to the density profile found in galaxies. 
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1 INTRODUCTION 

A galactic wind of dark matter particles (DMP) (see e.g. 
iBertone et al\ (2005)) flics through the Solar System and a 
part of it becomes captured due rotation of planets around 
the Sun. The capture process, dominated by Jupiter, is 
related to the nontrivial aspects of the restricted thre e- 
body problem (see e.g. IValtonen and Karttunenl |2005h ). 
We demonstrate that this process is des cribed by a sim - 
ple dynamical symplectic map ( see e.g. IChirikov! l|l979l ); 



iLichtenberg and Liebermanl (|l992l )) which allows to perform 
extensive numerical simulations of DMP capture. Our stud- 
ies show that the capture cross-section is much larger than 
the area of Jupiter orbit being diverging as an inverse square 
of DPM velocity in agreement with rec e nt an alytical esti- 
mates by iKhriplovich and Shepelvanskvl (|2009l ) . 

The dynamical map analysis allows to simulate DMP 
capture and ejection on the whole life time scale of the 
Solar System for 10 14 DMP being more efficient tha n 
the direct simulations of DPM dynamics by iPeteJ (2009). 
Our approach provides a DMP density distribution in 
the Solar System with other features of dynamics at 
present time after 4.5 billion years evolution of the So- 
lar System. This DM P distributio n is s imilar to those 
found in galaxies by iRuben et al\ |l980). The dynam- 
ics of DMP is shown to be chaotic having certain sim- 
ilarities with a chaotic comet motion in the Solar Sys- 
tem discussed bv lPetroskvl |l986l ); IChirikov and Vecheslavovl 
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J 19891); [Duncan et al\ J2005I); iDvorak and Kribbel l|l990j ); 
iMalvshkin and Tremaind (1 19991) . 

Following Bertone et al\ (|2005l ) we assume that in 
a vicinity of the Solar System (SS) the velocity dis- 
tribution of DMP has a Maxwell form f(v)dv = 
^54/ttu 2 /^ 3 exp(-3u 2 /2w 2 )d v with the average module ve- 
locity u « 220fcm/s. During a scattering of DMP on 
the Sun its rescaled total energy w — — 2E/mdV p is 
changed due to planetary rotatio n. The main contribution is 
given by Jupiter, as discussed bv lChirikov and Vecheslavovl 
i| 19891 ), and hence we base our studies on the case of one 
planet measuring DMP parameters in units of planet ra- 
dius r p and velocity v p taken as unity, DMP mas s mj = 



1. The studies of comet dyna m ics by Petroskvl (1 1986); 
IChirikov and Vecheslavovl l|l989h : iDuncan et al\ (|2005l ) in 
SS with one rotating planet show that it is well described 
by a symplectic map and thus a DMP dynamics over an 
extended orbit is also described by that type of map. 



2 DARK MAP DESCRIPTION 

This dark map has a form simi lar to the Halley map (see 
IChirikov and Vecheslavovl l|l989l )): 



-l =w n + F(x n ) , 



-l = x n + w. 



-3/2 



(1) 



where x n = t n /T p (modi) is given by time t„ taken at the 
moment of DMP n — th passage through perihelion, T p is 
the planet period. The second equation in ([T]) follows from 
the Kepler law for the DMP orbital period. The amplitude 
of kick function F(x) is proportional to the ratio of planet 
mass m p to the Sun mass Ms (F ~ m p /M s ) (see iPetroskvl 
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Figure 1. Le/f panel: dependence of kick function F(x) on 
Jupiter phase x for DMP orbit parameters shown by pluses: a) 
q = 0.11, 6 = 2.83, ip = 1.95 of the Halley comet case; here crosses 
show da ta for the Halley comet w i th all SS planets taken from 
Fig.l of IChirikov and Vecheslavovl Jl989l) ; b) q = 1.5, 9 = 0.7, 
ip = 7r/2; c) g = 0.5, 6 = 7r/2, ip = 0; curves show fit functions 
of numerical data marked by pluses. Right panel: dependence of 
maximal amplitude F max on q for a, b, c cases of left panel. 



l|l986h : IChirikov and Vecheslavovl (|l989l )). Its shape depends 
on DMP perihelion distance q, inclination angle 8 between 
the planetary plane (x,y) and DMP plane, and perihelion 
orientation angle ip. However, F(x) is independent of energy 
w for l/\w\ 2> r p — 1. Thus the dark map describes DMP 
dynamics for bounded and unbounded energies as well as its 
capture process corresponding to a transition from positive 
w < to negative energies w > 0. 

Our direct numerical simulations of Newton equations 
confirms this map description by F-function as it is shown 
in Fig. [1] for various val ues of q, 8, ip, including the Halle y 
comet case analyzed by IChirikov and Vecheslavovl (|l989h . 
In agreement with the theory of lPetroskvl ( 19861 ) the max- 
imum F m ax drops exponentially for q ^> r p so that only 
DMP with q < 2r p can be effectively captured. At q 3> 
r p we fi n d F ~ sin27ra in agreement with results of 
iPetroskvl (1 19861 ). The visible peaks in F max correspond to 
close encounters between DPM and planet happening at 
rather specific angles for q ^ r p . We will see later that 
such events give a small contribution in the capture cross- 
section a. In fact, F— function contribution comes from en- 
counter distances of the order of r p thus being much larger 
than t he radius o f the p lanet body r D . This analy ti cal re - 
sult of Petroskvl (Il986f); Chirikov and Vecheslavovl dl989h : 
IKhriplovich and Shepelvanskvl ( 20091 ); IShepelvanskvl (|2012l ) 
is in agr e emen t with the detailed numerical simulations 
by Peter J2009I) i nvalid a ting previous numeri c al stu dies of 
iGould and Aland (|200ll ); lLundberg and Edsicl (|2004D which 
considered contributions only from rb scale. 

Finally, we note that the dark map gives an efficient 
but approximate description. For the exact dynamics there 
is a slow variation of DMP o rbital momentum I and q = 
£ 2 /{2r p v 2 p ) and angles 9, ip fsee lDvorak and Kribbe] l |l99Ch ). 
However, the rate of these variations is rather slow being 
proportional to m p /Ms and does not affect significantly the 
chaotic diffusion i n energy. A lso numerical simulations of 
DMP dynamics bv lPeterl l|2009h point on a small global vari- 
ation of q. A similar situation appears in a microwave ion- 
ization of Rydberg atoms where it is known that the Kepler 
map in energy gives a good description of ionization process 
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Figure 2. Left panel: Dependence of capture cross-section a 
for Jupiter on DMP energy w; dashed lines show dependence 
a oc 1 / \w\; 1 / w 2 . Right panel: Dependence of rescaled captured 
number of DMP on energy w for Jupiter, Saturn and a model 
planet with m p /Mg = 0.004 (full, dashed, dot-dashed curves re- 
spectively) . 



of 3D-atoms as discussed by IShepelvanskvl (|2012l ). Also the 
DPM flow f(v)dv performs an averaging over all £, 0, ip val- 
ues and hence a variation of these parameters is averaged 
out. 



3 CAPTURE CROSS-SECTION 

In a scattering problem at infinity we have £ 2 = r|tjp|w 
with the impact scattering distance r 2 d — 2qr p /\w\. Hence, 
the capture cross-section at energy \w\ is a(w)/a p — 
(7r 2 r p |m|)~ 1 dd d<p J °° dqh(q, 6, ip), where h is a frac- 
tion of DMP captured after one map iteration from w < 
to w > 0, given by an interval length inside F(x) envelope 
at \w\ = const, a p = nr p . This fraction is determined from 
numerically computed F(x), as those shown in Fig. [1] via a 
continuous fit spline of function F(x). Using a grid with up 
to N g = 10 points in (q, 8, ip) volume we perform a Monte 
Carlo integration which gives cr(w) as a function of energy w 
for the case of Jupiter where the main contribution is given 
by |w| ~ w cap = m p v 2 /Ms « 10~ 3 . 

The dependence a(w)/(j p is shown in Fig. [5] For \w\ < 
Wcap we find a/a p w nMsw cav /m p \w\ in agreement with an - 
alytical estimates by IKhriplovich and Shepelvanskvl |2009), 
for \w\ > w cap we have a/a p ~ irMsw 2 ap /(m p w 2 ). The later 
regime describes contribution of close encounters which has 
a rapid decrease of a and hence gives a small contribution 
in the capture process. This conclusion is confirmed by the 
analysis of the differential number of captured DMP per 
time unit dN = a(w)n g v p f (w)d\w\ /2. Here n g is a galac- 
tic DMP density with a co rresponding mass de nsity p g — 
m d n g ~ 4-10" 25 g/cm 3 fsee lBertone et al\ (|2005l )) and f(w) 
is the velocity distribution function given above with \w\ — 
v 2 /v p at infinity. A number of DMP crossing the planet orbit 
area per time unit is N p = f 1 n g a p v p f(w)d\w\/2. 

The dependence of dN/N p dw on |tu| = v 2 /v 2 , pre- 
sented in Fig. [2] drops quadratically for \w\ > w cap show- 
ing that the contribution of close encounters is small. We 
note that dN/N p dw depends only on the ratio w/w cap that 
is confirmed by additional data obtained for Saturn and a 
model planet in Fig. [2] As a result the total number of cap- 
tured particles is N oc m v Ms i n agr eement with results of 
IKhriplovich and Shepelvanskvl (|2009l ) . 
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Figure 3. Left panel: The number N cap of captured DMP, as 
a function of time t in years, for energy range w > (dashed 
curve), w > 4 ■ 10~ 5 corresponding to half distance between Sun 
and Alpha Centauri System (black curve), w > 1/20 correspond- 
ing to r < 100 AU (red curve); DMP are injected at constant flow 
f(v) at all angles. Right panel: Top part shows density distribu- 
tion p(w) oc dN/dw in energy at time ts for DMP injection at 
all parameters q,8,ip (normalized as J pdw = 1), bottom part 
shows the Poincare section of the dark map for DMP injection at 
fixed parameters q, 9, ip of Fig. lb; inset shows density distribu- 
tion of captured DMP in w in log-log scale for parameters of Fig. 
lb (gray curve) and orbits of the main right panel injected at all 
parameters (black curve), a dashed line shows a slope -3/2. 



4 CHAOTIC DYNAMICS 

To determine the number of captured DMP N cap (t), in SS 
with Jupiter, as a function of time we model numerically 
a constant flow of scattered particles with energy distribu- 
tion dN s oc vf(v)dv per time unit. The injection, capture, 
evolution and escape of DMP is described by the dark map 
fl} with corresponding values of scattered parameters q,8,ip 
and corresponding to them F(x) function with the scatter- 
ing DMP distribution dN s oc dqdw (we use q < qmax = 4r p 
since above this value F max is very small). 

The scattering and evolution processes are followed dur- 
ing the whole life time of SS taken as ts = 4.5 • 10 9 'years. 
The total number of DMP, injected during time ts in the 
whole energy range ^ \w\ ^ oo, is N to t ~ 1-5 ■ 10 14 with 
Nh = 4 • 10 9 scattered DMP in the Halley comet range 
< \w\ sC w H « 0.005 (k = Ntot /Nh ~ 2u 2 /(3v%w H ) « 
3.8 • 10 4 , only DMP with \w\ < F max participate in dynam- 
ics). We used a random grid of initial q, 8, ip values with up 
to iVo = 4 ■ 10 J grid points and Ni injected orbits at each 
grid point with Nh = NoNi. 

The time dependence N cap (t) in Fig. [3] shows that ini- 
tially it grows linearly with time. This growth slows down 
after a time scale td ~ 10 7 years. For a finite SS region 
with w > 1/20 we see that there is a saturation of cap- 
tured number of DMP. Indeed, according to the results of 
IChirikov and Vecheslavovl (| 19891 ) for the Halley comet the 
scale td ~ 10 7 yr is a typical scale of diffusive escape of a 
comet or DMP from SS due to c haotic diffusion in energy. 
The a n alytical estimates given bvlChirikov and Vecheslavovl 
1 19891 ); IKhriplovich and Shepelvanskvl l^OoH )" also give a 
similar escape time. Thus after that time the injected flow 
is compensated by the escape process and we obtain a sys- 
tem in an equilibrium state with a fixed number of captured 
DMP with a certain energy distribution p(w). 

An example of such a distribution for a typical orbit 




Figure 4. Left panel: Radial density p(r) oc dN/dr at present 
time ts for SS with Jupiter averaged over all Nac DMP (the 
normalization is fixed as J Tp pdr = 1, r p = rj = 1). Right 
panel: Volume density p v = p/r 2 from the data of left panel, 
dashed line shows slope -2, horizontal line shows average density 
for r p /5 ^ r r p . 

parameters q, 8, tp at present time t = ts is shown in the 
inset of the right panel of Fig [3] (gray curve). There is a 
peak of density at small energies w < 0.002 where the orbital 
period is very long and chaotization is slow. For the range 
0.002 < w < w c h we have an approximate algebraic decay 
p ~ l/u> 3 / 2 which corresponds to the ergodic measure where 
DMP density is proportional to the orbital period T w ~ 
l/tii 3 / 2 . The chaotic diffusion to large energies is stopped by 
a critical invariant Kolmogorov-Arnold-Moser (KAM) curve 
which separates chaos region from integrable one at w = 
w ch . 

The analytical estimate of 

IKhriplovich and Shepelvanskvl (|2009| ) . based 



the Chirikov criterion (se e Chirikovl l|l979h : 
iLichtenberg and Liebermanl (|l992l )). gives for F max ~ 
5m p /Ms the value w c h ~ 0.3 that is in a good agreement 
with the case of Fig. [3] where w c h ~ 0.2. In a region 
w c h/2 < w < w c h we have stability islands, being well 
visible in the Poincare section, that gives significant fluctu- 
ations in density p(w). However, for w < w c h/2 the chaos 
component is homogeneous in the phase plane (w,x). This 
means that DMP are injected in the chaotic component of 
a chaotic layer around separatrix w = 0, and thus the DMP 
dynamics in SS is essentially chaotic. 



5 DENSITY AND MASS OF CAPTURED 
DARK MATTER 

To obtain DMP space density we consider Ntot scattered or- 
bits as described above. Their time evolution is described by 
the dark map |T} up to the present moment of time ts- We 
keep in memory the initial orbit parameters q,8,ip of cap- 
tured orbits. Then we consider only those with w > 4 ■ 10~ 5 
during the time interval 8ts/ts = ±10 -3 near time mo- 
ment ts collecting SNac ~ 6.2 ■ 10 6 orbits (while instan- 
taneously we have N A c ~ 3.3 ■ 10 5 ). For these 8N A c DMP 
their dynamics in real space is recomputed from their values 
of q, 8, ip, w, x during the time period of At ~ 100 Jupiter 
orbital periods using Newton equations. 

The radial density p(r) of DMP is obtained by averaging 
over 10 3 points randomly and homogeneously distributed 
over this time interval At for each of SNac orbits. The ob- 
tained normalized radial distribution p(r) in shown in Fig. [4] 
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Figure 5. Density of captured DMP in SS at present time t$- 
Top panels: DMP surface density p s oc dN/dzdr p shown at left 
in cross plane (0, y, z) perpendicular to Jupiter orbit (data are 
averaged over r p = \J x 2 + y 2 = const), at right in Jupiter plane 
(x,y, 0); only the range r ^ 6rj around the Sun is shown. Bot- 
tom panels: corresponding DMP volume density p v oc dN/dxdydz 
at left in plane (0, y, z), at right in Jupiter plane (x, y, 0); only the 
range |r| ^ 2rj around the Sun is shown. Color is proportional 
to density with yellow/black for maximum/zero density. 



with the corresponding average volume density p v = p/r 2 . 
It corresponds to a stationary equilibrium regime appearing 
at t 3> td when injection and escape flows compensate each 
other. The striking feature of the obtained result is that for 
r > r p we find p(r) ~ const, This means that the total DMP 
mass in a radius r grows linearly with r. 

According to the virial theorem such a profile gives 
a velocity of visible matter independent of radius v 2 n ~ 
Jo p{ r ')dr' /r ~ p(r), being similar to those found in galax- 
ies when th e DM P mass is dominant as discussed by 
iRuben et all j 19801 ) ; [Bertone et all |2005l ). Another impor- 
tant feature is that the DMP volume density p v remains ap- 
proximately constant for r < r p = r,j. However, for r > rj 
this density drops as inverse square distance from the Sun. 
Thus we find that a simple model of SS with one rotating 
planet is able to reproduce significant features of observed 
DMP density distribution in galaxies. 

Let us note that the radial density p(r) oc dN/dr is 
only approximately constant for r > r p . Indeed, a formal 
fit of data of Fig. [4] right panel in the range 2 < r/r p < 20 
gives p v ~ 1/r 13 with j3 = 1.53 ± 0.002. We can argue that 
this dependence can be understood from the ergodic mea- 
sure of effectively one-dimensional chaotic radial dynamics: 
dp ~ dN ~ pdr ~ J dtdw(dNfdw) ~ dt ~ dr/v r ~ y/rdr 
(assuming that dN/dw is peaked near w ~ as it is seen 
in the inset of Fig[3] and hence the radial velocity v r ~ 



if the DMP mass woul d be dominant, as i t is the case in 
galaxies as discussed bv lRuben etall Jl98CT). 

In fact the data presented bv lRuben et all (|l980t ) (see 
Fig.7 and Eqs.(l,2) there) are compatible with the depen- 
dence v m oc r 0,35 which is close to the above theoretical 
estimate. However, in SS the DMP mass is small compared 
to the visible matter and hence the case of galaxies should 
be analyzed in more detailed way using self-consistent condi- 
tions for DMP distribution which would modify the second 
equation in the dark map. Though the above arguments can 
be useful for analysis of DMP distribution at r 3> r p , in this 
work we perform the density analysis mainly inside the Nep- 
tune orbit where the radial density p(r) can be considered 
as approximately constant. 

From the data of Figs. [3] we determine the fraction 
Vac = N A c/Ntot ~ 2.2- 10~ 9 of DMP captured at time t s at 
energies w > 4- 10 and related fraction 7720 ~ 1-5- 10~ 41 at 
energies w > 1/20. From Fig. U we determine the fraction of 
Nac orbits in the volume r ^ 6r p with r) r a w 4.3- 10 -4 and in 
the volume r ^ r p with rj r i w 2.3- 10~ 5 . The DMP mass cor- 
responding to these fractions is obtained by multiplication of 
these fractions by the total mass of DMP flow passed in the 
corresponding range q ^ 4r p : M to t = J Q °° dvvf(v)ap g ts ~ 
69p g tskr p Ms /u « 0.9 • 10 _6 Ms where we use the cross- 
section a = 7rr| = 8ivkMsr p /v 2 for injected orbits with 
q ^ 4r p , k is the gravitational constant (u/v p w 17). Thus 
the mass of DMP with w > 4 ■ 10~ 5 is Mac ~ r)AcM to t ~ 
2 • 10- 15 Af s , and in a similar way the mass at w > 1/20 
is M 20 « rj2oM tot « 1.3 • 10" 17 M S . The mass Mac can 
be estimated as a mass of DPM with \w\ < wh absorbed 
by F ~ sin a; kick during the diffusion time td that gives 
Mac ~ v 2 w H tdMtot/(nu 2 ts) ~ W- 8 M tot ~ 10" 14 M S be- 
ing only by a factor 5 larger the above numerical value. 

The mass of DMP in the volume of Neptune orbit radius 
r < 6r p is M r6 = rj r6 M A c ~ 0.9- 10" 18 M S « 1.7- 10 15 3 and 
in the radius r < r p the DMP mass is M r i = r\ r \MAC ~ 4.6- 
Ms w 10 14 s. The average volume density of captured 
DMP inside the Jupiter orbit sphere r < r p — rj is pj — 
3M r i/(47rr 3 ) w 1.2 • 10~ 4 p 9 w 5 ■ 10- 29 g/cm 3 . Thus, the 
density of captured DMP is much smaller than the galactic 
DMP density. However, it is by a factor 4 ■ 10 3 larger than 
the equilibrium DMP galactic density p g u ~ 0.25p 9 /K 3//2 ~ 



1.4- 10" 



I cm taken in the energy range < \w\ < wh- 



1/y/r and p v ~ 1/: 



,3/2s 



Such a dependence 



would lead to velocity of visible matter v m oc ^fp oc 



.1/4 



The density distribution of captured DMP in SS is 
shown in Fig. [5] We see that the density decreases with 
r at r > rj in agreement with Fig. 3] A characteristic bulge 
is formed around the Jupiter orbit. A maximal local volume 
density is about 10 times larger than the average density pj 
inside r < rj. 



6 DISCUSSION 

For further studies it is desirable to take into account the 
contribution of other planets even if the results presented by 
IChirikov and Vecheslavovl (Il989l ) show that the main fea- 
tures of the dynamics are well described only by Jupiter 
contribution considered here. It is natural to expect that, 
as in the SS with one planet, the DMP dynamics in galax- 
ies is dominated by a few stars rotating around the cen- 
tral black hole and thus a constant radial DMP density 
behind Jupiter orbit found here should be typical for such 
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galaxies in agreeme nt with observational data discussed by 
iRuben et ai\ jl98tf ). 

In global our studies show that the average captured 
DMP density inside the Jupiter orbit is by a factor 10 4 
smaller than the galactic DMP density p g . The main rea- 
son for that is a small value of captured DMP energy 
w cap ~ mp/Ms ~ 0.005 which is very small compared to 
the dimensional galactic DMP velocity u/v p ~ 17. However, 
if we consider the galactic density in the capture energy 
range of < \w\ < w cap than we find that it is significantly 
enhanced by a factor 4 ■ 10 3 due to the capture process con- 
sidered here. Thus, the further analysis of chaotic capture 
process of dark matter in binary systems can bring interest- 
ing results. 

It would be also interesting to consider the inverse ion- 
ization process of DMP. According to the dark map |T|) the 
escape velocity square of DMP from a binary system of a 
star of mass m s rotating in a vicinity of a black hole of mass 
Mb is v\ ~ (m B /Mb)Vs. For a star moving in a vicinity of the 
Schwarzschild radius we may have the star velocity v„ ~ c/3 
and for the mass ratio m s /Mb ~ 0.01 we obtain the escape 
velocity of DMP Vd ~ c/30 w 10 4 km/s that is almost hun- 
dred times larger than the average galactic DMP velocity 
u ~ 200km/ s. Any other body of mass significantly smaller 
m 3 is ejected with a similar velocity that can generate com- 
pact wandering black holes crossing the universe at high 
velocity Thus, the stars on a distance of Schwarzschild 
radius from black holes can work as some kind of black hole 
accelerators generating high velocity DMP in the universe. 

We thank I.B.Khriplovich for useful discussions. A 
part of numerical computations has been performed at the 
mesocentre de calcul de Franche-Comte. 



REFERENCES 

Bertone G., Hooper D., Silk J., 2005, Phys. Rep. 405, 279 
Chirikov B.V., 1979, Phys. Rep. 52, 263 
Chirikov B.V., Vecheslavov V.V., 1989, Astron. Astrophys. 
221, 146 

Duncan M., Quinn T., Tremaine S., 1989, Icarus 82, 402 
Dvorak R., Kribbel J., 1990, Astron. Astrophys. 227, 264 
Gould A., Alam S.M.K., 2001, Astrophys. J. 549, 72 
Khriplovich LB., Shepelyansky D.L., 2009, Int. J. Mod. 

Phys. D 18, 1903 
Lichtenberg A. J., Lieberman M.A., 1992, Regular and 

chaotic dynamics, Springer, Berlin 
Lundberg J., Edsjo J., 2004, Phys. Rev. D 69, 123505 
Malyshkin L., Tremaine S., 1999, Icarus 141, 341 
Peter A.H.G., 2009, Phys. Rev. D 79, 103531; ibid. 79, 

103533 

Petrosky T.Y., 1986, Phys. Lett. A 117, 328 

Rubin V.C., Ford W.K. Jr. Thonnard N., 1980, Astrophys. 

Jour. 238, 471 
Shepelyansky D.L., 2012, Scholarpedia 7(1), 9795 
Valtonen M., Karttunen H., 2005, The Three-Body Prob- 
lem, Cambridge Univ. Press, Cambridge, UK 



© 2012 RAS, MNRAS 000, Offi] 



